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Abstract 

In a recent reformulation of three-dimensional new massive gravity (NMG), the field equations 
of the theory consist of a massive (tensorial) Klein-Gordon type equation with a curvature-squared 
source term and a constraint equation. Using this framework, we present all algebraic type D 
solutions of NMG with constant and nonconstant scalar curvatures. For constant scalar curvature, 
they include homogeneous anisotropic solutions which encompass both solutions originating from 
topologically massive gravity (TMG), Bianchi types //, VIII , IX, and those of non-TMG origin, 
Bianchi types VIq and VIIq. For a special relation between the cosmological and mass parameters, 
A = m 2 , they also include conformally flat solutions, and in particular those being locally isometric 
to the previously-known Kaluza-Klein type AdS-2 x S 1 or dS% x S 1 solutions. For nonconstant 
scalar curvature, all the solutions are conformally flat and exist only for A = m 2 . We find two 
general metrics which possess at least one Killing vector and comprise all such solutions. We also 
discuss some properties of these solutions, delineating among them black hole type solutions. 

PACS numbers: 04.60.Kz, 11.15.Wx 
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I. INTRODUCTION 



General relativity in three dimensions suffers from a number of barren properties: The 
usual Newtonian limit does not exist, there are no local dynamical degrees of freedom and 
no black hole solutions without a negative cosmological constant jl, 2] (see also j^] and ref- 
erences therein). With a negative cosmological constant, the theory admits black holes with 
an asymptotically anti-de Sitter (AdSs) behavior, which are known as Banados-Teitelboim- 
Zanelli (BTZ) black holes Q|. The BTZ black holes have played a profound role in under- 
standing the classical and quantum structures of the "simple" general relativity, albeit its 



jj. On 



par with this, intriguing devel- 



final quantum status continues to remain unclear 
opments have also been towards giving the theory dynamical degrees of freedom by adding 
higher-derivative terms to the Einstein-Hilbert (EH) action. Such a procedure generically 
introduces a topological mass parameter in the theory and entails the appearance of prop- 
agating modes with negative energy, known as ghosts. However, there are some fortunate 
exceptions in which cases the ghosts can be avoided. Among such exceptions, topologically 
massive gravity (TMG) occupies a central place jfj 7|. This theory is obtained by adding 
to the EH action a gravitational Chern-Simons term which is third-order in a derivative 
expansion. Nevertheless, the particle content of the theory is revealed in a usual way, by 
linearization of the field equations about a flat Minkowski spacetime. The resulting wave 
equation represents a single massive spin 2 mode of helicity either +2 or —2, depending 
on the sign of the topological mass parameter. This mode has positive energy only for the 
opposite sign (with respect to the conventional one in four dimensions) of the EH term in 



the total action. Thus, TMG is a unitary theory only 



and it is a parity- violating theory by its very nature |fj]. 

Recently, it was shown that there is a special case of TMG, determined by a critical 
value of the topological mass parameter, for which one can still allow the usual "right" sign 
of the EH term J In this case, the massive gravitons disappear and the theory admits 
BTZ black holes with positive mass. This in turn resolves the unitarity conflict between 
the bulk and boundary theories in the context of AdSs/CFT2 correspondence. That is, the 
theory of TMG in its special case becomes dual to a conformal field theory (CFT2) on the 
boundary (see Refs.jsl, 0] for details). This result has also renewed the interest in massive 
gravity theories in three dimensions on their own right and, as another fortunate exception, 



or the "wrong" sign of the EH term 



a new ghost-free theory of massive gravity appeared [10j. This theory is known as new 
massive gravity (NMG). It is a parity-preserving theory and is obtained by adding to the 
EH action a particular quadratic-curvature term. Though such an action gives rise to the 
field equations of fourth-order in derivatives, its canonical structure involves two propagating 
massive modes of helicities ±2. Moreover, the equivalence of this action, when expanding 
about the Minkowski vacuum, to the Fierz-Pauli action confirms unitary property of the 



ll|-ll3| 



theory. Subsequently, this property has also been confirmed in 

Further developments have been towards searching for exact solutions to the theory of 
NMG and studying its holographic properties. Shortly after the advent of NMG, it was 
found that the theory admits regular warped AdS3 black holes which are counterparts of 
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those appear in TMG Q. Continuing the hne of analogy wtth TMG S o,ution s 
the other constant scalar curvature solutions of NMG, such as AdS3-waves UM and Bianchi 



type IX homogeneous space solution [19[ were also found. A new class of solutions with 
nonconstant scalar curvature, which includes the black hole, gravitational soliton and kink 



type solutions, was discussed in 



can be found in 
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21] . Other studies of the issue of exact solutions to NMG 



22|-|26l|. The holographic properties of NMG exhibit some similarities with 



those of TMG: The bulk/boundary unitarity conflict still exists and its sensible resolution 
at some "chiral"_points turns out to depend crucially on the nature of asymptotic boundary 



conditions 



20 



27H29|. Therefore, finding pertinent stable vacua of the theory is of great 



importance and this motivates us to search for further exact solutions. 

In striving to undertake an exhaustive investigation of the exact solutions, the theory of 
NMG was recently reformulated in terms of a first-order differential operator (appearing in 
TMG and resembling a Dirac type operator) and the traceless Ricci tensor [30[. With this 
novel proposal, the field equations of NMG consist of a massive (tensorial) Klein-Gordon type 
equation with an "effective" curvature-squared source term and a constraint equation. This 
approach has a striking consequence for finding algebraic types D and N (in the Petrov-Segre 
classification) exact solutions to NMG. Remarkably, for these types of spacetimes, possessing 
constant scalar curvature, one can interpret TMG as the "square root" of NMG. This allows 
us to establish a simple framework for mapping of all types D and N solutions of TMG into 
NMG 30| . In other words, there is no need to solve in each case the field equations of NMG 
for a given metric ansatz as it becomes possible to generate solutions to NMG from the 
known solutions of TMG. Moreover, the novel proposal turns out to be a powerful tool for 



finding new types D and N solutions to NMG which do not have their counterparts in TMG. 
Some intriguing and 81m p.e examples of these "non-TMG origin" so.utions were given fl 
In the meantime, a universal metric for all type N solutions of NMG was found in [3 II ]. 

The purpose of this paper is to present all algebraic type D solutions of NMG with both 
constant and nonconstant curvatures. In Sec. II we begin by describing a theoretical frame- 
work which includes reformulation of NMG in terms of a first-order differential operator, 
appearing TMG and acting on the traceless Ricci tensor. In Sec. Ill we introduce an or- 
thonormal basis of three real vectors (one timelike and two spacelike vectors) and describe 
the properties of type D spacetimes in NMG, in terms of these vectors and their covariant 
derivatives. Here we also prove two defining mathematical statements for these spacetimes. 
In Sec. IV we present the exact solutions with constant scalar curvature, which include homo- 
geneous anisotropic space solutions of both TMG and non-TMG origins. Next, we present 
conformally flat solutions with constant scalar curvature, for which the cosmological and 
mass parameters are related by A = m 2 . In Sec.V we discuss the solutions with nonconstant 
scalar curvature and possessing at least one Killing vector. They all are conformally flat 
and require A = m 2 by their very existence. We find two general metrics which comprise all 
such solutions. We also discuss some properties of these solutions, among them are black 
hole type solutions. 



II. THEORETICAL FRAMEWORK 



For our purposes in the following it is convenient to begin with TMG. The theory is 



described by the action 



7] 



(1) 



where R is the three-dimensional Ricci scalar, A is the cosmological constant, /i is the mass 
parameter and the Chern-Simons term is given by 



Cos = \e» vp vi P (d^ + 1 r^) . 



(2) 



Here e^ vp is the Levi-Civita tensor and the quantities T^r are the usual Christoffel symbols. 
The field equations of TMG, obtained from this action, are given by 



Gnu + ^9au H C 



/!/' ) 



(3) 



where = R^ u — \ g^R is the Einstein tensor and C^ v is the Cotton tensor given by 

C^ = e^(R^-- A g v pR] (4) 
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and the semicolon stands for covariant differentiation. We note that this tensor is a sym- 
metric, traceless and covariantly conserved quantity 

It is not difficult to see that with the traceless Ricci tensor 



S^u — R^v — — g^uR , (5) 



equation (j3J) can be put in the form 

S,„ + -C, u = 0. (6) 
Next, we introduce a first-order differential operator If) which is defined as 

= I (V^c^ + e^^p) , (7) 
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where is a symmetric tensor. It is straightforward to show that with the subsidiary 
relation 

= (*„■%., (8) 

equation ([7]) takes the most simple form 

^<V = e M a/3 3W. (9) 
Since condition (jSJ) is fulfilled for the Schouten tensor 

®^ = R ^~ ^9»» R i ( 10 ) 

it is easy to verify that 

c^ = -#^ = -#<V. (ii) 

With these relations the field equations of TMG, as seen from equation ([6]), acquire the form 

IpS, v = iiS, u . (12) 

We note that this equation resembles the Dirac equation If)^ a = 1a_ B ^^b = for a 
massive spinor field. In the case under consideration, the operator Jp acts on the traceless 



Ricci tensor, IJ)S^ V = D^ L V pS n p, and therefore it looks like the "cousin" of the Dirac 



operator (see also Refs. [30l. l3l|). 

Using equation f fl2|) . it easy to show that the secondary action of the operator Jp yields 
the Klein-Gordon type equation 



(jp 2 - /i 2 ) <v = o 



(13) 



It is also important to note that with the differential operator Jp we have the following 
identity 



pa c fiv 



This identity can be verified by using the relation 

along with the fact that in three dimensions the Riemann tensor is given as 

(R 
Rp[a9j3]u — Ru[a9l3}p ~ 9p[a9/3]u 

where the square brackets stand for antisymmetrization over the indices enclosed. 
There also exists another important identity 



(14) 



(15) 



(16) 



{ip 2 s^ = c, v c» v - s^ip 2 s l 



pv 1 



(17) 



that is obtained by a straightforward calculation of the quantity S^Ip 2 S^, using the defi- 
nition in ([7]), with equations ([9]) and ( TTTT) in mind, as well as the identity in ( fl4l) . 
We now proceed to NMG whose action has the form 1(| 



(18) 



where A is the cosmological parameter, m is the mass parameter and the scalar K is given 
by 
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K = R, V W V --R 2 . 

pf „ 



(19) 



In a recent work 



301 ] . it was shown that the field equations, which follow from this action, 



can be written in terms of the square of the operator Jp , resulting in the Klein-Gordon type 
equation 



(^ 2 -m 2 )^ = T, 



p,u 1 



(20) 



where T^ v can be thought of as the effective source tensor. We have 

R 1 

T^ u = S^pS p u — — S^v — — g^ySapS . (21) 

Equation (120]) is also accompanied by the subsidiary (constraint) equation 

S^S^ + m 2 R - ^ = 6m 2 A . (22) 

It is remarkable that such a description of NMG greatly simplifies the search for exact 
solutions to the theory. It is instructive to demonstrate this using the case of maximally 
symmetric solutions to TMG, for which the Cotton tensor in ([6j) vanishes identically and we 
have S^ IU = 0. Then equation ( 120]) is trivially satisfied, whereas equation ( 122]) yields 



A = 2m 2 (1 ± yjl - X/m 2 j , (23) 

where we have used the fact that R = 6A. Thus, with the cosmological constant adjusted 
according to this relation, every maximally symmetric solution of TMG (including its iden- 



tification with a BTZ black hole [4|) can be mapped into two inequivalent solutions of NMG 

HQ. 

The most striking feature of this description emerges for algebraic types D and N space- 
times. It is straightforward to verify that for these spacetimes the following relation 

T^ u = kS^, (24) 

holds, where k is a function of the scalar curvature. With this relation, equation ( 120]) reduces 
to the form 

(IP 2 -^)S, U = 0, (25) 

where 

/j, 2 = m 2 + K. (26) 

We see that for type N spacetimes and type D spacetimes with constant scalar curvature 
this equation is equivalent to that given ( TTB"]) . Thus, comparing equations ( 1T2"]) . ( TT3"]) and 
( 125]) we conclude that for these types of spacetimes TMG can be interpreted as the square 
root of NMG. This is a remarkable fact and it paves the way for mapping all types D and 
N exact solutions of TMG into NMG by means of an algebraic procedure that adjusts the 



physical parameters of the corresponding solutions in both theories 



30| . On the other hand, 



it is clear that equation ( 125|) admits types D and N exact solutions of non-TMG origin as 
well. (The square root does not always exist!) Some simple examples of such solutions were 
given in jjjf]]. Meanwhile, in 31 1 we managed to find a universal metric which comprises all 



type N solutions of NMG. Here we wish to continue the line of works 
an exhaustive investigation of all type D solutions of NMG. 



30] and [31], giving 



III. TYPE D SPACETIMES 



Let us introduce an orthonormal basis of three real vectors {t^, s M , v^}, satisfying the 
relations 



s„.s h 



(27) 



with all other products vanishing, such that the spacetime metric can be written in the form 



9nv 



(28) 



In three dimensions, depending on whether the one-dimensional eigenspace of S^, is timelike 
or s pac elike, one can distinguish types D t and D s spacetimes, respectively [32| (see also 



33 



34]). For type D t spacetime, the canonical form of the traceless Ricci tensor is given by 

S^u = P (fiW + 3Vf) , (29) 



whereas, for type D s spacetime we have 



(30) 



where p is a scalar function. In what follows, we shall perform a detailed analysis for type 
D t spacetimes. The same analysis can be performed for type D s spacetimes as well, either 
directly or simply by an appropriate analytical continuation of the type D t results. 

Clearly, using expression ( )29|) in equation ( 122|) one can express the function p in terms of 
the scalar curvature, the mass and cosmological parameters, m and A. Thus, we have 

R 2 



6p = 6m X — m R + 



24 



(31) 



On the other hand, taking into account expression (|2"9"|) in equation (I2~4"j) and using the result 
in equation (|26|) . we find that 

R 



jj, = m — p 
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(32) 



We turn now to the Cotton tensor in (j3j). It is straightforward to verify that the most 
general form of this tensor, written in terms of the orthonormal basis vectors, must have the 
following representation 

- Vfj,v u ) + 2cv( p s v ) + 2ht^s u ) , (33) 

where the round brackets denote symmetrization of the indices enclosed. The coefficients a, 
b, c and h are scalar functions and we have gauged out the term proportional to t^v u ) using 
the invariance of equations (|28|) and (|29|) with respect to rotations in the (s M , f At )-plane. 
Using this representation in equation ( 1T4"|) . we find that 

i^S^Y = 3phv». (34) 

Next, taking the divergence of equation ( )25|) and comparing the result with equation (|34|) . 
we arrive at the relation 

Sj{^)-p + \ \i 2 R, a = 3phv a , (35) 
o 

where we have used the contracted Bianchi identity 

S^ V = \R^. (36) 

Further it is convenient to consider in ( 132]) the cases of fj, 2 ^ and /i 2 = separately: First, 
we begin with the case fi 2 ^ 0. Contracting both sides of equation (135]) with s M , we find 
that 

s^R 4l = 0. (37) 

That is, the scalar curvature is constant along the basis vector s M . Meanwhile, the contrac- 
tion of equation f )35|) with yields 

h = -^—v^R.., (38) 

where we have used equations (I3TI) and (132|) . We note that the same contraction procedure 
with gives a trivial result. Next, we need to find the representation for the covariant 
derivative of the vector t M . For this purpose, we substitute expression (1291) into equation 
( ITTl) and taking into account the relations 

tp £[ii/pS ^ i &fj, £pvpt Up £p,vpS (39) 



together with the Bianchi identity in (136]) . we compare the result with the representation of 

the Cotton tensor given in (l33l . After straightforward calculations, we find that 

a 2b c g , / 

W = g V % - — s (m ^ } + — (fys„ - + — (s M s„ + + — v^t v , (40) 

where 



9=Y^=-{p' + h)^ t (41) 

/= (p'-^)^» (42) 

the prime stands for the derivative with respect to the scalar curvature. In obtaining ( 14 ip 
we have also used expression (J29l) in equation (136]) with the subsequent contraction of the 
result with 

We are now ready to prove that the following statement: Suppose that the scalar function 
h = 0. Then, for [i 2 ^ 0, all type D solutions of NMG must possess constant scalar curvature. 

The proof will be given by contradiction, assuming that R ^ const. First, we note that 
for h = equation (|38|) yields 

v»R ;t t = . (43) 
Moreover, comparing equations (j!4p and fl34l) we obtain that 

C p[tl S/ = 0. (44) 

This relation is trivially satisfied in the TMG case, whereas it is a governing geometrical 
relation for the NMG solutions. It is important to note that with h = in ( 1331) . one can 
always gauge out either b or c due to the rotational symmetry of the spacetime metric (1281) 
in the (s M , w M )-plane. Therefore, in the following we can set c equal to zero, without loss of 
generality. 

We also note that with relations (1371) and (|4"3"|) . one can always write down the simple 
decomposition 

R;» = Xt», (45) 

where x is some function and \ ^ as R ^ const by our assumption. Using the fact that 
= 0? we take the covariant derivative of this equation. As a consequence, we have 

Xtfcu] +t\nXv] = 0. (46) 
10 



Substituting now f|40l) into this equation and contracting the result with v^s v ', we find that 
a = 0, in addition to c = 0. Thus, equation fHOj) reduces to the form 



2b g . , . 



Meanwhile, the Cotton tensor in (1331) becomes as 

C^v = b (s^ - v^v u ) . (48) 

The value of b in these equations can be fixed by using ( ITTj) along with ( 1251) and ( 134|) . We 
find that 

b 2 = 3pV- (49) 

Next, we calculate the divergence of expression (l4~8l) . Taking into account the property 
= 0, and contracting the result with the vectors s M and t> M , we have 

V M = ° ; ^ ;M = 0, (50) 



where we have also used relations (1371) and (1431) . 

With equations ( 1471) and ( |50|) . as well as using the properties of the orthonormal basis 
vectors given in (|27j) . it is straightforward to write down the explicit expressions for the 
covariant derivatives of the basis vectors s M and t> M . They are given by 

b g 

= -—t^v u + —t^s u + (v^U , (51) 
3p 3p 

v ix;v = — 2^ t^ s u + ^ ^^i/ — C s ^u > (52) 

where £ is a scalar function. 

Finally, the use of these expressions in the Ricci identity 

x ^\v il — {%^). v = + x>1 1 (53) 

written for the vectors s M and v^, gives us ( = as well as the relation 

(^ 2 )'(l + ^)^ ;At = 0. (54) 

We recall that fi 2 ^ 0. Furthermore, using equations (15Tj) and (1521) . it is easy to show that 
the expression in the round brackets does not vanish identically. Thus, from equation ( |54|) 
(see also (157|) and ( 115]) ) it follows that the scalar curvature i? is constant. This contradicts 

11 



to our initial assumption that R 7^ const and completes the proof of the statement made 
above. 

A similar analysis shows that this statement holds for type D s spacetimes as well. How- 
ever, in both cases an interesting question is that what happens if one drops the condition 
h = 0. It turns out that in this case one can still prove the same statement, provided that 
the spacetime admits a hypersurface orthogonal Killing vector. On the other hand, the 
general case with h ^ becomes very involved and unclear. 

We turn now to the case /1 2 = 0, which is equivalent to 

P = ™ 2 - Y^, (55) 

and taking this into account in equation fl3T|) . we find that 

A = m 2 . (56) 

Meanwhile, from the field equation (125]) it follows that Jp 2 S^ v = 0. Using this in equation 
(GUI, we see that h = 0. With these results, from equation ( FT71) it follows that 

C ia/ C ta ' = 0. (57) 

Substituting ( 1331) into this equation, we find the relation 

3pV + b 2 = , (58) 

which implies that both a = and 6 = 0. Therefore, 

= , (59) 

as can be seen from equation (We recall that c is gauged out as h = 0) . Thus, with 

a 2 = 0, or equivalently A = m 2 , all type D solutions of NMG must be conformally flat. 



IV. SOLUTIONS WITH CONSTANT SCALAR CURVATURE 

We shall now discuss all type D solutions with constant scalar curvature, focusing first 
on the case /i 2 7^ 0. These are homogeneous anisotropic solutions which can be formally 
divided into two categories: Solutions of TMG origin and solutions of non-TMG origin i.e. 
those which do not have their counterparts in TMG. 

12 



We begin by noting that for R = const, the functions h, g and /, as seen from equations 
( 138|) . (j41j) and fH2l) . vanish. Then, equation (T4T)]) can be written in the form 

a 2b 

W = 3V V^W' ( 60 ) 
Similarly, equation ( 1331) reduces to the form 

Cfu, = aSf,,, + b (s^s v - VpVy) . (61) 

Contracting now equation ([25]) with tH v and using equations (jUj) and (fTTj) . after some alge- 
braic manipulations, we obtain that 

(2p a t p;CT + C/V) = 2p^ 2 . (62) 

Substitution expressions (!60|) and (1611) into this equation yields the relation 

a 2 + |j=/. 2 . (63) 

On the other hand, writing down the Ricci identity (|53|) for the vector and contracting 
the result with t u , we find the relation 

a 2 -- 2 =9p-*-R. (64) 
2 

Since p is constant, as can be seen from (EH), comparing equations ( 1631) and ( 1641) it is easy 
to see that a and 6 are constants as well. Furthermore, with equation ( IBTj) it is not difficult 
to see that for b ^ 0, we again arrive at the relations given in (!50l) . 



A. Solutions of TMG origin 

When 6 = 0, the Cotton tensor in (IBTj) becomes as 

C^y = aS^y , (65) 

with 

a = ±/i . (66) 

Comparing this equation with those given in (fTTl) and f[T2"j) . and taking into account (1131) and 
(1251) . we deduce that all type D solutions of TMG can be mapped into NMG. The resulting 
solutions, as can be seen from equations (|32|) and (|64l) . are characterized by 

m 2 17 . . 

13 



It is straightforward to show that the use of this value of p in (131]) and (132|) results in the 
adjusting relations, earlier obtained in {30], for mapping of the TMG solutions into NMG. 
For these solutions equation ( 1501) takes the form 

W =3 e F %=3 - > ( 68 ) 
where in the last step we have used the relations in ( 139|) . That is, the vector tp is the Killing 



vector (see also 



13, 



I32|). 

We now want to give an elegant classification of these solutions by a dimensional reduction 
on this Killing vector, where the sign of the scalar curvature of the two-dimensional subspace 
(the factor space) plays a crucial role. We define the projection operator onto the factor 
space perpendicular to tp. It is given by 

h\ = 5\ + tH v , (69) 

and Wyf = 0. We also define the derivative operator Dp with respect to the two-dimensional 
"spatial" metric hp U , and the associated Riemann tensor r fJl u \ T . We have 

D u Vp = h\h\ V x , a , r\ Xr Vp = 2D [T D X] V U . (70) 

With these definitions, the relation between the spacetime Riemann tensor and the Riemann 



tensor of the factor space is derived in a usual way (see for instance, Ref. 35[). After 
straightforward calculations, we have 

r fii/al3 = h u h P a H p (R\ T p a + 2 t T - t \ tp- a + tp- : \ t r;cr — t p;T t\- a ) . (71) 

Taking the projection of equations (fTBj) and f )68|) onto the factor space and using the result 
in (I71j) . we find that 

r^vafi = f 2p + — - — J (h m h v p - hupKa) ■ (72) 
From this expression it follows that the scalar curvature of the factor space is given by 

r=?ifl-^, (73) 
20 5 ' K ' 

where we have also used equations (152]) and (!67|) . 

Next, we need to find the covariant derivatives of the basis vectors Sp and Vp. Using 
equation f )68|) . together with the properties of the orthonormal basis in (j2TJ), it is not difficult 

14 



to show that the most general expressions for these derivatives are given by 

s^-u = - - tpV,, + AvpU + BvpSv + Cv^v v , (74) 

v n;u = g t^i s u ~ As^ty — Bs^s v — CS^Vy , (75) 

where the functions A, B and C need to be determined. Projecting these expressions onto 
the factor space, by means of (1691) . and comparing the result, for certainty, with that obtained 
for the metric of two-dimensional hyperbolic space 

ds 2 2 = — (dQ 2 + sinh 2 9 d<p 2 ) , k 2 = - T - , (76) 
k ^ ' 2 

we find that 

5 = 0, C=kcoth6. (77) 
Substitution of these quantities in equations (1741 and ( 1751) yields 

s^ = C, v M = 0- (78) 



On the other hand, using the Ricci identity (15311 for the vector t> M , with equation ( 17811 in 
mind, and contracting the result with v u , we find that A = —a/3. Thus, equations (174|) and 
f5j) take their final forms given by 

s^-u = ~ t| (t^Vp + Vffiv) + Cv^Su , (79) 



v 



)1\V 



- (tpS,, + SpU) - Cs^Sy . (80) 

It is not difficult to show that the associated Lie brackets are given by 

[M] = o, M = 0, [v i8 ] = ~t + Cv. (81) 

It is convenient now to choose a coordinate system in which the Killing vector is given as 
t = d T . Then, the use of (176]) and (IHTI) enables us to specify the remaining vectors in the 
form 

k / 2a \ 

s=kde > v= ^{ d *-w coshe v- (82) 

Meanwhile, the corresponding dual one-forms result in the spacetime metric 

2a \ 2 1 

— cosh^j +- 



/ 2a \ 2 1 

ds 2 = - ( dr + — cosh 6 d<p\ + — (d6 2 + sinh 2 6 1 



(83) 
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We see that this metric has an apparent singularity at k = 0. However, by redefining the 
coordinates as 



r ^ r "f^' 9 ^ k9 > ( 84 ) 



one can put it in the form 



2 



r 2 ( r 4a . . 2 k9 , \ 
as — — I dr + — ^ smh — dtp I 

sinh 2 k9 , l9 

+ d9 2 + — p — d0 2 , (85) 

which is suitable for taking the limit k — > 0. 

It is also worth noting that in terms of the left invariant one-forms of SU(1,1), parame- 
terized by the Euler angles, 

(J\ — — sin ip d6 + cos xjj sinh 9 d<f) , 
(i2 = cos ip d0 + sin ■?/> sinh d0 , 

o"3 = + cosh 6 d(j) , (86) 
and with r = (2a /3k 2 ) ip, the metric in ( |83l) takes its most simple form given by 

^--(wf^ + h^+rf- (87) 

In this case, instead of (!8~T]) . we have the following Lie algebra 

[s,t) = -—v, [v,t] = —s, [v,s] = -—t. (88) 

We recall that the quantities a and k 2 are as given in ( )66|) and ( |76|) and they are determined 
by equations ( |73l) and ( IBTl) together with (|3TT) and (1321) . We see that, in general, depending 
on the sign of the scalar curvature r of the two-dimensional factor space, we have three 
possible spatial geometries. Namely, a sphere for r > 0, a hyperboloid for r < and a 
flat space for r = 0. Accordingly, the associated spacetimes are homogeneous anisotropic 
solutions: (i) the solution of Bianchi type IX, or with SU(2) symmetry, (ii) the solution 
of Bianchi type VIII, or with SU(1, 1) symmetry, (iii) the solution of Bianchi type II, or 
with the symmetry of the Heisenberg group. 

We note that the type D s counterpart of spacetime (1H71) can be obtained by taking t — >■ is 
and s — > it in the above expressions (see also Eqs. (l29|) and (!30|) ). In this case, Bianchi type 
// remains unchanged, whereas Bianchi types VIII and IX goes over into each other. 
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B. Solutions of non-TMG origin 



We turn now to the case of b ^ in equations (1501) and (1BT|) . This results in solutions 
which are only inherent in NMG. Clearly, in this case the vector t M is no longer the Killing 
vector, but it is still divergence- free. Furthermore, as we have mentioned above, the vectors 
s M and are divergence-free as well. With these in mind, using equation (160]) and the 
properties of the orthonormal basis vectors given in (1271) . we find that 



v K v = ^ (a ~ H tfiSu - ipSfJv , (90) 



where ip is a scalar function. Using these relations in the Ricci identity ( 1531) for the vectors 
s M and v^, and taking into account equation (150]) . we find that 

^ = 0, P=-f- (91) 

Next, we calculate the action of the operator Jp on the Cotton tensor in ( IST]) . by means of 
equations ( TTTj) . ( |89j) and ( 1901 . Then, comparing the result with equation ( l25j) . we find that 
a = 0, whereas 6 is given by the same relation as in (l4"9]) . (See also Eq. (l63~]) ). 
On the other hand, from equations (!3T|) . (l32j) and (164"}) it follows that 

4 , m 2 , , 

p = ~i5 ' T" (92) 

Finally, using these results in equations ( 160|) . (189|) and (1901) it is straightforward to show 
that we have the following Lie algebra 

[s,t] = ^2/5mv, [v,t] = ^2/5ms, [v,s] = 0. (93) 

It is convenient to choose a coordinate system in which the Lie algebra admits the following 
representation for the basis vectors 

t = d T , s = ^ (e~^ mT d x + eV^ mT d y ^j , 

v = -j= [e-^ mT d x - eV^^dy^ . (94) 

Then, the spacetime metric can easily be written down by using the associated dual basis. 
We have 



= _ dr 2 + e 2^2/ 5 mr dx 2 + e -2^2/ 5 mr dy 2_ (g5) 
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This solution, as can be seen from the Lie algebra in (193]) . is of a homogeneous anisotropic 
spacetime of Bianchi type VIq, or with E(l, 1) symmetry. 

Again, the corresponding type D s solution can be obtained by making the replacement 
t — > is and s — > it. In this case, we have the Lie algebra 



[t,s] = -fi/Smv, \p,s] = j2/hmt, [v,t]=0, (96) 
which enables us to specify the basis vectors in the form 



s = d x , , v = sin(v/2/5 mx) d T + cos(J2/5mx) d y 



t= -cos(^2/5mx)d T + sm(\/2/5mx)d y . (97) 

The associated dual basis results in the solution given by 

ds 2 = cos (2 ^2/5 ms) (-dt 2 + dy 2 ) + dx 2 + 2 3^(2^2/5 mx) dtdy . (98) 

This is a homogeneous anisotropic spacetime of Bianchi type VIIq, or with E(2) symmetry. 

To complete this subsection, we consider now type D solutions with constant scalar 
curvature, which correspond to the case n 2 = 0. For these solutions, we have the special 
relation between the cosmological and mass parameters, as given in (156j) . and the Cotton 
tensor in (I6T!) vanishes as both a and b are zero. Using these facts in equation (160]) . we see 
that 

t KV = . (99) 

It follows that t = d T is a hypersurface orthogonal Killing vector of constant length. With 
this in mind, we employ the Ricci identity ( 1531) for the vector t^ and obtain that 



P=f- (100) 
Using this value of p in equations ( 1551) and ( 1721 . we find that 

R = r = 4m 2 . (101) 



Clearly, type D t spacetime metric can now be written as 

1 

k 



ds 2 = -dr 2 + — Ud 2 + sinh 2 9 d<p 2 ) , (102) 



where 

k 2 = -2m 2 = -2A . (103) 
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Meanwhile, type D s solution is obtained from this metric by performing the coordinate 
changes r — > i<p and <fi — > ir. We have 



ds 2 



? + — (c^ 2 -sinh 2 #dT 2 
k 1 v 



(104) 



It is straightforward to show that this metric is locally isometric to the Kahrza-Klein solution 
AdS 2 x S 1 found earlier in [14J (see also Ref. 20]). We note that the spacetime metric f 1 1 2 j) 



has also a "cousin" with spherical spatial section (the flat case is trivial). One can show 
;hat the type D s counterpart of the latter is locally isometric to the dS 2 x S 1 solution of 



201 ] . We recall that similar solutions with a hypersurface orthogonal Killing vector absent 



in TMG due to the no-go theorem of 



V. SOLUTIONS WITH NONCONSTANT SCALAR CURVATURE 

Now we discuss all type D solutions which possess nonconstant scalar curvature. We 
recall that these are conformally flat solutions for which C^ v = 0(a = 6 = c = 0). Using the 
traceless Ricci tensor ( 129|) in the Bianchi identity (|36|) and contracting the result with the 
vector t M , we obtain that 

V M = °> ( 105 ) 

which in turn implies that g = 0. In obtaining these results we have also used equations 
()4TT) and (l55j) . With these in mind, from equation ( 1401 we have 

V = — V- ( 106 ) 
v 

In what follows, we will restrict ourselves to the case when the spacetime under consideration 
admits at least one Killing vector £. For further convenience, we use the invariance of metric 
(1281) with respect to rotations in the (s M , t> M )-plane and represent the Killing vector in the 
form 

6* = oAp. + , (107) 

where a and (3 are scalar functions to be specified later. Then, due to the same rotations 
we have 

t KV = — — V- + s„t v , (108) 
p V 
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instead of (I106p . From the Bianchi identity ( 136)) it follows that 

eS,v'" = 0. (109) 
Substituting equation f[29~j) in this expression and using the relation (11051) . we find that 

0. (110) 

On the other hand, the use of the fact that 

£^ = 0, (111) 
where £^ is the Lie derivative along the Killing vector, gives us the relations 

(ap). M s" = 0, (ap) J X = 0- ( 112 ) 
From equations f II 1 j) and (11121) it follows that 

a = aop' 1 , (113) 

where ao is a constant. 

Using the properties of the orthonormal basis in (127)) along with expression (11081) and the 
Killing equation 

=0, (114) 



we obtain the relations 



and 



P P 

v mv = tpt v H r-^ s^s v , (116) 

p p 



= s m Am = °- ( 117 ) 



Next, we need to obtain the determining equations for unknown functions /3 and p. 
They are obtained by employing the Ricci identity ([5*5)1 for the basis vectors s M and v M , 
respectively, as well as using equations (j!08j) . dl 15[) and (I116p . After some manipulations, 
we arrive at the "oscillator" type equation for /3, 

v»v»/3. KV -k 2 /3 = 0, (118) 
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as well as at the following set of equations for p . 



s^(p-%., u = 3 + k V 1 - ^ v»(p-% , (119) 



v*v v {p- x ) wv = S + tfp- X , (120) 



v^ip-%.^0. (121) 



We note that k 2 is the same as that given in f ll03p . It is straightforward to show that the 
Killing vector in f llOTj) commutes with the basis vector v in the sense of their Lie bracket. 



That is, we have 

[£,«] = 0, (122) 

which can easily be verified by means of equations f ll08j) . f II 1 5 j) and (I116p . This fact allows 
us to choose the vectors £ and v as 

£ = d z v = d x (123) 

in a coordinate system (x ,y ,z). With this in mind, we turn to equation ( 11 18j) whose general 
solution is given by 

/3 = h{y) sinh(A;x) + / 2 (y) cosh(A;x) , (124) 

where f\{y) and f2{y) are arbitrary functions. Alternatively, the use of the reparametrization 
invariance of the vector fields £ and v with respect to the coordinate freedom x — > x + g\ (y) 
enables us to reduce this solution into three different forms given by 

P = u(y)smh(kx) , (125) 
(3 = u(y)cosh(kx) , (126) 

/3 = u(y)e kx , (127) 

where the function u>(y) will be specified below. 

We turn now to equations (11071) and fll 13j) and consider the cases a = and a ^ 
separately. 
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A. The case ao = 



In this case, the Killing vector, as seen from equation (I107p . becomes proportional to the 
spacelike basis vector s. That is, we have 

s = ^d*- (128) 

Next, using the relations 



[£,*] = 0, [ V ,t] = ^t, (129) 

where the first relation follows from equation (jllip and the second one is obtained by means 
of equations ( I108p and ( JUS]) , we fix the timelike basis vector as 

t = pd y . (130) 

With this in mind, from the first relation of (JTTTJ) , we find that u = const in equations 
( I125j) - (|127j) . Using now the dual one-forms to the basis vectors in equations (11231) . d 128[) 
and (I130p we arrive at the spacetime metric 

ds 2 = -p~ 2 dr 2 + dx 2 + /3 2 dy 2 , (131) 

where we have relabeled the coordinates as y = r and z = y. We see that this metric in 
general is characterized by two functions (3 = (3(x) and p = p(r,x). The first function is 
given by solutions (I125p - (I127I) . whereas the second function is determined by the system of 
differential equations fl 1 19j) -( TT2Tj) . which admits the general solution 

p-^-^ + ftrW, (132) 

where /(r) is a smooth function. Performing an analytical continuation of solution f II 3 1 j) 
by t — > iy and y — > ir, we arrive at its type D s counterpart. It is straightforward to show 
that the resulting spacetime metrics for the explicit forms of j3 given in (H25p - fT27]) . after 
appropriate coordinate changes, represent the black hole, gravitational soliton and kink type 
solutions with one Killing vector, respectively. Meanwhile, for / = const the second Killing 
vector appears as well, and these solutions reduce to those with two commuting Killing 
vectors, found earlier in 2l| . Below, we shall focus only on the black hole type solution with 



one Killing vector. 
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Let us take = smh(kx) in metric (11311) . Then, passing to the coordinates 3/k 2 — r = 
Mcosh(ya;), where r is a radial coordinate and M is a constant, r — » i<ft and y — » ikMr, 
we arrive at the metric 

(ir 2 

o?s 2 = — fc 2 (r — r + ) (r — r_) c?r 2 + 



k 2 (r — r+) (r — r_ 

2 7 i 2 



+ [r- (r + + r_)/2 + F(^)f^, (133) 

where the radii of outer (r + ) and inner (r_) horizons are given by 

r± = -^±M, (134) 

and -F(</>) is an arbitrary function. This metric can be interpreted as describing a "gen- 
eralized" black hole type solution, which is asymptotically AdS spacetime with A = 2m 2 
(m 2 < 0). For r + = r_ = r , it corresponds to an extremal black hole with one Killing 
vector, which was earlier described in [30(. It is also easy to see that for F((f>) = const, 
this solution goes over into a standard new black hole metric (with two commuting Killing 
vectors) of [2o|,|21|. 

In light of this, it is also interesting to ask whether the solution in fl 1 3 1 [) admits another 
limit with two Killing vectors. It turns out that the answer is affirmative, but the Killing 
vectors are no longer commuting. In fact, solving the equation rju.^ = for the putative 
Killing vector n, we find that 

V = d T - ^ d x + b yd y , (135) 

provided that the metric functions are given by 

/ = e 6or , f3 = e kx . (136) 

Here bo is a constant of integration. We see that the Killing vectors n and £ do not commute 
for 6 7^ 0. Clearly, similar analysis is also true for the type D s counterpart of (11311) . 



B. The case ao ^ 

We begin with calculating, in addition to (I122p . the other Lie brackets of the vectors £, 
s and v. Using equations (I107p and ( II 15)) and taking into account equation (11171) it is easy 
to show that 

[£,«] = 0, (137) 
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whereas, the use of equations (I115P and f 1 11 6 j) yields 



[s,v] = -^8. (138) 

With these Lie brackets and with equation (11231) . one can show that the vector s admits the 
simple representation 

s = ±d y , (139) 

and lo = const in equations (H25p -f lT2T|) . In obtaining these results we have also used the 
first equation in (11 ITj) and the coordinate freedom z ^ z + g 2 (y) . 

Next, using the above representation of the vectors £ and s in equation (I107p . we see that 
the timelike basis vector t is given by 

t = P (d Z -dy), (140) 

where we have set = 1, for certainty. It is easy to see that the associated dual one-forms 
enable us to write down the spacetime metric in the form 

ds 2 = -p~ 2 dr 2 + dx 2 + p 2 {dy + udrf , (141) 

where we have passed to the new coordinates p —> (3/uj , r — >■ z and y — > ycu. We recall 
that p = p(x) is determined by the solutions in (11251) - 11271) . and p = p(x, y) is given by the 
solutions of equations (I119I) - (I121I) . Thus, we have the following three pairs of solutions 

p' 1 = -— + c cosh(kx) + {c x e iky + c 2 e~ iky ^ sinh(A;x), 

p = smh(kx), (142) 

p^ 1 = —— + c smh(kx) + (cie ky + c 2 e~ ky ^ cosh(kx), 

p = cosh(A;x), (143) 



p- 1 = - p + c e kx + ci (k 2 y 2 e kx + e~ kx ) + c 2 ye kx , 



3 

p = e kx , (144) 



where Co , c\ and c 2 are constants of integration. The type D s counterpart of solution (j!4ip 
is obtained by making the coordinate changes r — >■ iy and y — >■ ir. 

24 



It is important to note that metric (11411) with one Killing vector is intrinsically different 
from that in ( JTHTj) even for u = 0. In particular, this can be seen from the fact that for 
metric (1131)) . for which a = 0, the Killing vector in (I1U7P is an eigenvector of (|2T?j) . i.e. 
Sf^iu — V^v However, this is not the case for metric ( 1141 j) . where ^ and we have 
S^tiv — V^v~ 3ao^- On the other hand, for the vanishing constants C\ and c<i an additional 
Killing isometry appears in (jl4ip . and d y becomes the second Killing vector, commuting 
with the first one, £ = d T . With these two commuting Killing vectors, solutions (1 1 3 1 p and 
(I14ip become locally isometric to each other. 

In analogy with the case of (1131)) . it is not difficult to show that for the special forms of 
the metric functions, solution (1141)) still admits a second Killing vector. This vector does 
not commute with the Killing vector £. Again, solving the associated Killing equation, we 
find the explicit form of the Killing vector. It is given by 

n = e k(ulT+y) [tanh(fcx) d y - d x ] , (145) 

provided that the metric functions are as in equation (11431) with cq = and c<i = 0. We 
note that for u = 0, this vector commutes with the Killing vector £. Then, it is not difficult 
to show that the resulting metric with these commuting Killing vectors is locally isometric 
to that in (11311) with / = const and with (3 = e kx . It is also worth to note that the Killing 
vector, similar to that in (I145p . does also exist for the solutions in (11421) . again with the 
vanishing cq and C2- On the other hand, there is no such a Killing vector for the solutions in 
(11441) . Finally, it is certainly interesting to know what is the relation between metrics (1131)) 
and (1141)) in their special limits with the Killing vectors in (1135)) and ( 11451) . respectively. 
This issue requires further investigation. 

In summary, we note that the spacetime metrics in ( 1131)) and (11411) . both admitting one 
Killing vector, represent all type D solutions of NMG with nonconstant scalar curvature. 
We recall once again that all these solutions are conformally flat and they do exist provided 
that the special relation A = m? holds. 



VI. CONCLUSION 



This paper completes our programme of the exhaustive investigation of types D and N 



exact solutions of NMG, which was first begun in 



30] and then continued in 3JJ. Here 
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we have presented an exhaustive set of solutions, which includes all type D solutions with 
both constant and nonconstant curvatures. As in the previous cases, this was achieved in 
the framework of a novel proposal that amounts to reformulation of NMG in such a way 
that the field equations underlying the theory acquire a remarkably simple form. Namely, 
they reduce to a massive (tensorial) Klein-Gordon type equation which is accompanied by 
a constraint equation as well. The "mass term" in this equation plays a distinguished role 
in many aspects of the delineation of exact solutions. In particular, for its constant value 
being achieved only for type N and type D (with constant scalar curvature) spacetimes, 
there exists an intimate relation between the theories of NMG and TMG; the latter can be 
understood as the square root of the former, thereby opening up the way for mapping all 
the associated solutions of TMG into NMG. 

We have proved that for the nonzero value of the mass term, all type D solutions of NMG 
must have constant scalar curvature, whereas for the vanishing mass term all the solutions 
turn out to be conformally flat, possessing both constant and nonconstant scalar curvatures. 
Introducing an orthonormal basis of three real vectors and using the gauge freedoms pro- 
vided by the associated Lorentz symmetries, we have given a detailed description of type D 
spacetimes in NMG, in terms of these basis vectors and their covariant derivatives. 

With type D solutions of constant scalar curvature, we have shown that all the solutions 
fall into two classes: (i) homogeneous anisotropic solutions which consist of those having their 
counterparts in TMG and those of being only inherent in NMG, (ii) solutions with vanishing 
Cotton tensor, i.e. conformally flat solutions. Using the dimensional reduction procedure 
on a Killing vector, we have given an elegant classification of the homogeneous anisotropic 
solutions of TMG origin in terms of the scalar curvature of two-dimensional subspace. We 
have obtained the most compact expression for the spacetime metric and established a 
universal Lie algebra for the associated basis vectors. Depending on the value (positive, 
negative or zero) of the two-dimensional scalar curvature, appearing in the Lie algebra, there 
exist three possible spatial geometries and the associated spacetimes are of Bianchi types 
IX, VIII and II, respectively. We have also obtained all homogeneous anisotropic solutions 
of non-TMG origin, which are of Bianchi types VIq and VIIq spacetimes. It is important 
to note that these solutions require a special relation between the cosmological and mass 
parameters (A = m 2 /5), unlike the Bianchi types of TMG origin. As for conformally flat 
solutions, they exist only for A = m 2 and possess a hypersurface orthogonal Killing vector of 
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constant length. We have discussed these solutions as well, emphasizing that some of them 
are locally isometric to the previously-known Kaluza-Klein type AdS2 x S 1 or (IS2 x S 1 
solutions. 

With type D solutions of nonconstant scalar curvature, we have found two general metrics 
which admit at least one Killing vector and comprise the entire set of such solutions. All 
these solutions are conformally flat and require A = m 2 by their very existence. We have 
discussed the special limits of these metrics with two commuting Killing vectors. In the 
latter case, the resulting metrics become locally isometric to each other, recovering all the 
previously-known solutions in the literature. We have shown that there also exists a special 
limit of both general metrics when two noncommuting Killing vectors appear. The question 
of whether in this case the resulting metrics are locally isometric or not is unclear and it 
requires further investigation. As an illustrative example, we have briefly discussed global 
properties of the solutions with at least one hypersurface orthogonal Killing vector, focusing 
on black hole type solutions. 

In this paper, we have given all the solutions in their local form. It is certainly of 
great interest to perform the complete global analysis of all type D solutions presented here, 
especially in the case with nonconstant scalar curvature. It is also of great interest to extend 
the results of this paper to other theories of 3D massive gravity. Using our approach, this 
can be done for "general massive gravity", for an extension of NMG by adding a parity 
violating Chern-Simon term, as well as for the theory of NMG with higher-order curvature 
invariants. We hope to return to these issues in our future works. 
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